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Abstract 

We present a stochastic computational model to study the mechanism of signalling be- 
tween a source and a target ionic transporter, both localized on the plasma membrane (PM) 
and in intracellular nanometre-scale subplasmalemmal signalling compartments comprising 
the PM, the sarcoplasmic reticulum (SR), Ca^"*" and Na"*" transporters, and the interven- 
ing cytosol. We refer to these compartments, sometimes called junctions, as cytoplasmic 
nanospaces or nanodomains. In the chain of events leading to Ca^"*" influx for SR reloading 
during asynchronous Ca^^ waves in vascular smooth muscle (VSM), the physical and func- 
tional link between non-selective cation channels (NSCC) and Na''"/Ca^''' exchangers (NCX) 
needs to be elucidated in view of two recent findings: the identification of the transient 
receptor potential canonical channel 6 (TRPC6) as a crucial NSCC in VSM cells and the 
observation of localized cytosolic [Na"*"] transients following purinergic stimulation of these 
cells. Having previously helped clarify the Ca^"*" signalling step between NCX and SERCA 
behind SR Ca^^ refilling, this quantitative approach now allows us to model the upstream 
linkage of NSCC and NCX. We have implemented a random walk (RW) Monte Carlo (MC) 
model with simulations mimicking a Na^ diffusion process originating at the NSCC within 
PM-SR junctions. The model calculates the average [Na"*"] in the nanospace and also pro- 
duces [Na"*"] profiles as a function of distance from the Na^ source. Our results highlight 
the necessity of a strategic juxtaposition of the relevant signalling channels as well as other 
physical structures within the nanospaces to permit adequate [Na"*"] build-up to provoke 
NCX reversal and Ca^"*" influx to refill the SR. 

Keywords: calcium oscillations, calcium signaling, sodium transient, vascular smooth mus- 
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1 Introduction 



We are not used to thinking of ionic sodium (Na"*") as a signalling ion, despite its known 
relevance for vascular disease. On the other hand, the importance of the second messenger 
Ca^^ in signalling cell function is undisputed. There is however recent experimental evidence 
supporting the idea that this important species has, at least in vascular smooth muscle 
cells (VSMC), an almost equally important signalling partner in Na"*"!!], [2], particularly 
in signalling events that are allowed by the juxtaposition of signalplexes and transporter 
carrying membranes at nanometric distances from one another. This observation, and others 
outlined below, as well as a wealth of accumulated knowledge on Na^ and Ca^^ transporters 
in VSMC has prompted us to take a more quantitative look at the question of Na^-related 
signalling in PM-SR cytoplasmic nanospaces (also referred to as nanodomains or junctions): 
nanometre-scale signalling compartments comprising the PM, the sarcoplasmic reticulum 
(SR), Ca2+, Na+ and other ionic transporters (channels, exchangers and pumps) and relevant 
signalplexes therein, and the intervening cytosol. 

In this article we present a quantitative model aimed at elucidating the mechanism of 
selective (or site-specific) signalling between a source ionic transporter and a target ionic 
transporter, both of which are localized on the same membrane and are part of a nan- 
odomain. We developed this model using a stochastic method based on the simulation of 
ionic diffusion by random walks (RW) within nanospaces modeled according to experimental 
observations. Here we concentrate on the specific example of a Na^ transporter, typically a 
NSCC and a NCX as its target. Generally, during these events Na"*" entry via a NSCC gen- 
erates a large [Na^] gradient, which, in turn, enables reversal of a NCX (presumably in the 
vicinity of the NSCC) and consequent NCX-mediated Ca^"*" entry into the subplasmalemmal 
nanodomain. A few examples highlighting the biological importance of this intra-membrane 
system and its link to pathogenic mechanisms are (table [1] summarizes these cases): 1) 
A critical subplasmalemmal step in the Ca^^ signalling cascade giving rise to VSM cell 
contraction following G-protein coupled receptor adrenergic stimulation [3] — where blocking 
NCX reversal causes attenuation and elimination of Ca^^ oscillations due to impairment of 
SR refilling; replacement of [Ca^^ji oscillations with a tonic Ca^"^ signal causes a dramatic 
decrease in smooth muscle force development |1] ; in pulmonary artery SM, upregulation of 
the NCX and Ca^+ entry via reverse NCX action is considered one of the mechanisms behind 
elevated [Ca^"*"] in idiopathic pulmonary arterial hypertension patients|5]; 2) Receptor acti- 
vated Na^ entry may also induce NCX reversal in endothelial cells (EC) and this, in turn, 
gives rise to selective Ca^'''-stimulated eNOS activity and NO production; eNOS derived NO 
is an important physiological vasodilator agent and is accepted as an independent marker 
of vascular health[6l [^; in EC NCX operating in forward mode is also important in regu- 
lating both [Ca^"*"]! and [Ca^^]sR; this suggests that local [Na''']i, besides the cell membrane 
potential and the equilibrium potential of the NCX, have a role in the regulation of forward 
NCX too[8l[9]; 3) In nerve terminals, a transient [Na^] increase linked to a tetanic pulse of 
the action potential can reverse the NCX to induce presynaptic potentiation; this observa- 
tion implies a role for the NCX in synaptic facilitation and has consequences for short term 
memory from reverse NCX malf unct ion [TOj [2]; 4) In skeletal muscle NCX plays an impor- 
tant role in Ca^^ homeostasis, operating in forward as well as reverse mode in that function; 
Na^/Ca^^ exchange is involved in the control of muscle fatigue and there are reports sup- 
porting the notion that the beneficial role of external Ca^"*" in protecting slow-contracting 
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blood flow 
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blood flow 
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rev and fwd/eNOS activity 


NO production regulation 
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nerve 


rev/tetanic pulse 


short memory 






of action potential 






skeletal muscle 


rev/Ca^"*" influx 


muscle fatigue 


m 



Table 1: Sampling of systems in which NCX modes are linked to pathologies. Legend: 
rev=reverse mode (Ca^^ influx); fwd=forward mode (Ca^^ efflux); ACaW=Asynchronous 
Ca2+ waves; PASM=pulmonary artery smooth muscle; EC=endothelial cells. 



soleus muscle against high-frequency fatigue depends mostly on Ca influx through reversal 
of the NCX[Tl]; 5) although the role of NCX in heart is still poorly understood, both its 
Ca^^-efflux and influx modes are observed in cardiac myocytes and and the latter is likely a 
consequence of subplasmalemmal Na"*" elevations [12]. 

All of the above emphasize that modulation of Ca^"*" signaling via Na"*" and Na'''/Ca^''' 
exchange is of great clinical relevance in areas such as hypertension, chronic heart failure 
and possibly cerebral and skeletal muscle malfunction. 

This study sheds new light on a few new scientiflcally interesting key factors regarding the 
"workings" of intracellular signalling nanospaces. We investigate the role, and importance, of 
having a conflning surface, namely another lipid membrane, facing the membrane where the 
source and target transporters belong. This emerges as an important feature of nanodomains, 
as previous observations had hinted [H] . but it would appear that it alone cannot be 
entirely responsible for the generation of a sufficiently large local [Na^] transient elevations. 
We find that to understand how sufficient [Na''"] can be built up within the nanodomain to 
trigger NCX reversal and Ca^^ signalling downstream, it is important to account for other 
factors such as the role of physical obstructions to Na"*" diffusion and the possible organization 
of these obstructions. 



2 Methods 

2.1 Electron microscopy 

Details of the electron microscopy have been described previously [15]. The primary fixative 
solution contained 1.5% glutaraldehyde, 1.5% paraformaldehyde and 1% tannic acid in 0.1 M 
sodium cacodylate buffer that was pre-warmed to the same temperature as the experimental 
buffer solution (37 °C). The rings of rabbit IVC were fixed at 37 °C for 10 minutes, then 
dissected into small blocks, approximately 1 mmxO.5 mmxO.2 mm in size and put in the 
same fixative for 2 h at 4 °C on a shaker. The blocks were then washed three times in 0.1 M 
sodium cacodylate (30 min). In the process of secondary fixation, the blocks were post-fixed 
with 1% Os04 in 0.1 M sodium cacodylate buffer for 1 h followed by three washes with 
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Figure 1: To-scale model nanospace used in the simulations. The separation between the 
two surfaces (cyan and blue) is 20 nm. 

distilled water (30 min). The blocks were then further treated with 1% uranyl acetate for 1 
h (en bloc staining) followed by three washes with distilled water. Increasing concentrations 
of ethanol (25, 50, 70, 80, 90 and 95%) were used (10 min each) in the process of gradient 
dehydration. 100% ethanol and propylene oxide were used (three 10 min washes each) for the 
final process of dehydration. The blocks were infiltrated in the resin (TAAB 812) and then 
embedded in molds and polymerized in an oven at 60 °C for 8 h. The embedded blocks were 
sectioned on a microtome using a diamond knife at the thickness of 80 nm. The sections were 
then stained with 1% uranyl acetate and Reynolds lead citrate for 4 and 3 min, respectively. 
Images were obtained with a Phillips 300 electron microscope at 80 kV. 

2.2 Simulations 

Simulations of Na"^ diffusive motion from a source transporter are based on an implementa- 
tion of a Monte Carlo random walk (RW), along the lines of methodology previously employed 
for Ca^"*" transport simulations |14j. The model nanospace in which the simulations of Na"*" 
diffusion take place arises from observed physical features and properties of the essential 
elements for nanodomain signalling. These are either obtained from our laboratory's studies 
or from the available literature [16] and they are essentially the TRPC6 cytosolic 'radius' 
and FmaxIHl, the typical dimensions of intracellular nanospace ultrastructure, estimates on 
the number of pillars, Na"*" diffusivity in cytosol, expected [Na^] necessary for NCX reversal 
(see section 13.11) and during localized Na"*" transients (LNats) observed in [1] . Fig. [1] is a 
to-scale representation of the geometry of the model nanospace used in the simulations. In 
the simulations, particles representing Na"*" performs a RW on a cubic lattice with spacing 
s = 0.2 nm; initially picked as an approximation to the expected Na"*" mean free path in 
water (in turn, as an estimate to the mean free path in cytosol), we also carried out tests for 
effects of varying this parameter in a previous article and revealed no substantial influence 
on the results[Tl]. The RW time step r was chosen by running several simulations letting 
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particles cover a predetermined straight line distance and recording the number of RW 
steps taken to cover d. From diffusion theory, the total time t taken by a random walker 
in three dimensions to cover the distance d is t = d'^ /{GDj^eas), where -Dmeas is the measured 
diffusivity of Na"*" in muscle cytosol[IH]- The quantity t/N was our choice for r and its value 
is 10-11 s. 

The boundary conditions in our simulations are as follows. At the PM and SR membranes 
we implemented reflecting conditions: ions arriving at one of those surfaces during their RW 
are reflected back into the nanospace. At the edges of the model nanospace conditions 
are perfectly absorbing: ions reaching the lateral boundary of the junction are considered 
absorbed by the cytosol external to the nanospace, lost from the junctional population and 
no longer contributing to the PM-SR nanospace concentration, [Na''']ns- As explained later, 
we also positioned a number of obstacles (we refer to them as pillars) to Na"*" motion spanning 
the distance between the membrane in the nanospace. Pillars behave like elastic scatterers 
for Na+. 

Simulation code is written and tested in the C programming language on a computer 
running a Linux operating system. After testing and troubleshooting, programs are recom- 
piled and run in one of the WestGrid computing nodes|19]. The pseudo-random number 
generator we used is the algorithm gsl_rng_ml9937 of the GNU Scientific Library [201 El], 
since it has sufficient randomness and quality requirements for our purposes. The plots with 
simulation results were produced using the "freely distributed plotting utility" gnuplot[22]. 

3 Results 

3.1 Model foundation 

This laboratory has previously experimentally established that a [Na^] transient from a 
NSCC enabled the generation within a PM-SR junctional nanospace of a sufficiently high 
[Na"*"] gradient to cause NCX reversal[3]. A quantitative estimate of the level of such burst 
can be obtained by comparing the equilibrium potential of the NCX, -Encx, with the mem- 
brane potential, Vu, using typical values for vascular smooth muscle cells. These quantities 
are linked by the equation -Encx = Vm (where -Encx = 3i?Na — 2i?ca) between the electro- 
chemical potential of Na"*" (-^Na), of Ca^"*" (-Eca) and the membrane potential (Vm)[9]. (In 
general, Ec^+ = i?T/(zF) ln([C^+]out/[C'''''~]in), where R is the universal gas constant, T the 
absolute temperature, C^"*" is a 2;-valent cation, and F is Faraday's constant.) NCX reversal 
occurs when -Encx < ^m- We can study this inequality for both resting and activated cell 
conditions by a plot like the one in Fig. [2j Using for an activated cell = ~20 mV, 
[Na+]o = 140 mM, [Ca^+Ji = lO/xm, [Ca^% = 2 mM, R = 8.3 J/(mol K), T = 310 K, and 
F = 9.65 X 10^ J/(V mol), we observe that, during activation, a [Na"*"] transient of the order 
of 30 mM or greater is necessary to cause NCX reversal. In this exercise, we have used an 
estimated value for [Ca^'*"]! = 10/im, as that is approximately the lowest value [Ca^~i']i we 
expect from observations cited earlier [T^. 

Now, equipped with (a) the fundamental observation of localized [Na"*"] elevation tran- 
sients in full agreement with the values suggested by the study of Fig. I2(lj . (b) better 
knowledge of the identity of NSCC as TRPC6[Il [23], (c) the basic idea that the presence of 
intracellular nanospaces is necessary for this signalplex to be complete, we propose a model 
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Figure 2: Equilibrium NCX potential in the case of a resting (red curve) or activated (green) 
smooth muscle cell compared with the corresponding membrane potential, Vm- 

to investigate the role of strategic placement of transporters with respect to each other, as 
well as of a confining membrane and other ionic diffusion limiting structures in the generation 
within these nanospaces of sufficiently high [Na^jns to permit NCX reversal. 

The model nanospace used for the study is illustrated in Fig. [H The dimensions expressed 
therein are based on high quality EM images showing that the PM-SR separation in these 
nanospaces is remarkably uniform and about 20 nm. Lateral extension of these closely 
apposed PM-SR regions is approximately 400 nm (see [H]). 

3.2 Random walk simulations: bare PM-SR nanospace 

The simplest model nanodomain we studied consists of a shallow cylinder-shaped volume of 
height h = 20 nm and radius R = 200 nm, as in Fig. [1] Na"*" entering the nanodomain via an 
NSCC are represented as particles doing a RW based on a given diffusivity D = 600/im^/s, 
corresponding to that of Na"*" in muscle cytoplasm as reported in [18]. For simplicity, there is 
only one Na"*" source positioned at the centre of the PM side of the nanospace (in section H] we 
will elaborate further on the issue of the number of Na"*" sources) . The simulation programs 
output the computed [Na''']ns rise above resting level as a function of time, thereby giving 
an average [Na"*"] increase in the nanospace. (From now on [Na''']ns denotes the increase in 
cytoplasmic Na"*" concentration in the nanodomain between PM and SR.) Results from this 
set of simulations are shown in Fig. [3l left panel. This plot helps establish the time scale 
after which we can consider that the [Na"*"] has reached a steady-state level. We can observe 
from the graph that after approximately 100 us the concentration has reached a plateau after 
having increased from zero during an initial transient, as we have explained in the previous 
section. Having established a time scale for the formation of the approximately maximum 
level of [Na''']ns, we compute and study the concentration profile inside the junction, by 
plotting [Na''']ns as a function of the distance from the Na"*" source. The graph in the right 
panel of Fig. [3] illustrates a representative result. 
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Figure 3: Left panel: [Na''']ns vs t for simplest form of nanospace as in Fig. [U right panel: 
[Na''']ns vs distance from Na''" source placed at the centre of the PM side of the nanospace. 



3.3 Random walk simulations: randomly distributed obstacles in 
the nanospace 

Evidently, this simplest incarnation of the model is inadequate to describe the generation 
of Na^ transients of the observed size [I], since at steady state elevation of [Na''']ns hovers 
around 2 x 10~^ M or about three orders of magnitude less than the observed values of 15-20 
mM. We need to consider other nanospace features emerging from our ultrastructural images 
that may be responsible for a larger increase in the [Na''~]ns- Barring artificially changing the 
value of Na"*" diffusivity, [Na''']ns can be "forced" to increase if the ions were able to dwell 
longer in the nanospace than they are in the simple version of the system analyzed so far. 

There is convincing evidence suggesting the existence of structures spanning the width of 
the nanospace and which could constitute an impediment to the free diffusion of Na^ [2| [2l] . 
Our own observations confirm the existence of electron opaque "pillars" in transmission 
electron microscopy images like the one in Fig. HI The size and abundance of these electron 
opaque structures compares well with the electron dense "bridges" observed by Devine and 
collaborators in the early '70s|24j. 

Keeping the overall geometry of the model nanospace the same (Fig. [1]), we have therefore 
implemented a number of junction spanning structures in the form of cylindrical pillars, 
having estimated their size from several images like the one in Fig. HI From those same 
images it is possible to approximate the percentage junctional volume occupied by those 
structures and, in turn, an approximate number of them expected per nanospace. In a 
series of simulations, we have represented up to 200 pillars randomly distributed within 
the nanospace (Fig. [5]), and then simulated the Na"*" diffusion by a random walk within the 
pillar-populated nanospace. In this case too, we let the simulations run for a time sufficiently 
long to ensure that a steady-state level for the [Na^jns was established. Results are reported 
in Fig. El In all simulations involving random positioning of pillars, to minimize bias from 
the particular random pillar distribution, we have take average values of the computed [Na^] 
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Figure 4: Transmission electron microscopy image showing nanospace bridging electron 
opaque structures (white circles). Tissue: rabbit inferior vena cava. SR=sarcoplasmic retic- 
ulum; PM=plasma membrane; ECS=extracellular space. 

over 10 different random pillar distributions. These are the values plotted in the graphs we 
present in this article. 

The results of this series of simulations indicate that having some form of impediment to 
ionic motion in the junction does produce the effect of increasing the [Na^jns and of changing 
its profile to one that decays more slowly with distance from the Na"*" source. We ensured 
that this effect was not simply a consequence of the decreased nanospace volume due to the 
presence of the pillars by plotting the steady-state [Na"*"] computed with different numbers 
of pillars in the junction (blue dots in Fig. [7]) and comparing it with an increase in [Na^] 
merely due to reducing the nanospace volume by the volume of the pillars (red line in Fig. 
[7j). The plot in Fig. [7] demonstrates that ion collisions with pillars do indeed have a role in 
forcing Na''' to dwell longer in the nanospace. 

3.4 Random walk simulations: non-randomly distributed obsta- 
cles 

Clearly, the presence of obstacles to diffusion has an effect of increasing [Na''']ns, however 
the values we obtain this way are not yet comparable with those measured during the local 
[Na"*"] elevation transients. Other junctional features need to be accounted for in order to 
understand the mechanism giving rise to such high Na"*" transients necessary to reverse the 
NCX and observed by Poburko and coworkers[T]. We considered the hypothesis that nature 
might place these obstacles "strategically" rather than randomly, in a neighbourhood of a 
Na"*" source, so as to favour the generation of the gradients needed to drive the signalling 
chain. The rationale is that while a random set of pillars does show an ability to retain ions 
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Figure 5: To-scale model nanospace including about 200 10-nm radius randomly placed 
cylindrical pillars spanning the distance between the two surfaces. 

in the nanospace longer and therefore allow higher concentration build up, it does not do 
it efficiently enough to quantitatively account for the observed [Na''']ns transients. We then 
ran some simulations in which pillars are placed in a circle around a Na"*" source in such a 
way that we can control the porosity of this pillar fence to the passage of random walking 
ions. (Imagining to stand where the Na"*" source is, the circle of pillars would appear as a 
20-nm-high set of slabs surrounding the source itself, with thin rectangular gaps between 
the slabs; with this in mind, porosity is defined as the ratio of the entire surface area of the 
gaps to that of the slabs plus gaps.) The configuration of our model nanospace in this case 
is shown in Fig. [HI as a two dimensional x,?/-projection. 

Representative results from these simulations are reported in Fig. [9l The dramatic 
increase in the [Na''']ns in the vicinity of the source caused by this type of obstacle configu- 
ration is immediately evident. [Na''']ns is in this case of the same order of magnitude as the 
observed localized [Na^] transient elevation phenomenon [1]. Note furthermore that in this 
configuration the model suggests that the time for the [Na''']ns to reach steady-state is much 
longer than in the case of random distribution of pillars. The more scattered data in the 
plot of the right panel in Fig. [9] are due to the fact that the system in this case has yet to 
reach equilibrium. 

4 Discussion 

The stochastic computational model presented herein attempts to give a quantitative descrip- 
tion of the mechanism behind the observed localized [Na"*"] transients observed in VSMC[T]. 
Based as much as feasible on experimental observations of the physical and physiological 
features of the intracellular nanospaces in which these transients are hypothesized to occur, 
the model results lead us to conclude that [Na"*"] can build to sufficiently high values in 
PM-SR junctions to give rise to the observed transients. 

Three main steps lead to the fundamental hypothesis behind this work. The first stems 
from earlier work by this laboratory elucidating the sequence leading to asynchronous Ca^'*' 
oscillations in VSMC[T1[T3]. Succinctly, a large external Na"*" infiux causes the reversal of the 
NCX and consequent Ca^"*" entry to refill the SR after SR-Ca^"*" release via IP3R channels 
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Figure 6: [Na^jns within a model nanospace containing a number of pillars occupying about 
30% of the volume, as in Fig. [51 Left panel: computed [Na"'"]ns as a function of time; 
the values of [Na^] in this graph are an average across the entire nanospace. Right panel: 
concentration profile within nanospace. The red curve and dots are the same data displayed 
in Fig. [31 



upon cell stimulation provokes contraction. The second is the observation of two main 
features of the localized [Na^] elevation transients [ij: they appear as a punctate pattern 
on the periphery of VSMC (with puncta having a given time course), and their peak [Na^] 
values are comparable to the [Na"*"]! values necessary to cause NCX reversal (see Fig. [21 
and relevant text). Thirdly, our previously published model supports the idea that NCX 
reversal-mediated Ca^^ entry in PM-SR nanospaces introduces sufficient Ca^^ to refill the 
SR during asynchronous Ca^"*" waves [T^. Based on these three points, the hypothesis we set 
out to study with our model is that the observed localized [Na"*"] transients occur in PM-SR 
nanodomains (or nanospaces or junctions), in other words, that in those junctions, due to 
Na^ infiux via a TRPC6 channel, [Na"*"] can attain levels, that cause Ca^^ entry via NCX 
reversal. 

Simulation results from our simplest version of model (Fig. [3]), namely, a PM-SR 
nanospace filled with cytosol (represented in our simulations by the diffusivity of Na"*") with 
only one TRPC6 channel as Na"*" source suggest that this simple view is not adequate to de- 
scribe the formation of [Na"*"] transients. This is mainly due to the large value of the diffusion 
coefficient of Na"*", which is about three times that of free Ca^"*" in cytosol, and the fact that 
there is no observed buffering effect of Na^ in cytosol[18]. Recent [Na"*"] measurements in iso- 
lated rabbit ventricular myocytes by Despa and Bers ([251 [26]) also suggest that the effective 
diffusivity of Na"*" may be much slower than the one found by Kushmerick and Podolsky|18j. 
although no explanation as to the mechanism behind it was suggested. In [25] , among other 
things, Despa and Bers measured endogenous Na^ buffering in cytosol and found it negligible 
when compared to that of other important species like Ca^^, for example. It is well known 
that Ca^"^ buffering lowers its diffusivity dramatically [271 [28]. thus possibly contributing to 
easier generation of Ca^"^ gradients in confined spaces like the PM-SR nanospaces. Presum- 
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Figure 7: Comparison of steady-state Na"*" calculated only accounting for the volume effect 
of pillars in the nanospace (red line) and computed from the simulations (blue dots). 



ing the slight Na buffering effect observed in cardiac cell cytosol is mirrored in smooth 
muscle cells, it clearly could not provide sufficient slowing down of Na"*" diffusivity to aid 
local transient creation. It seems instead ever more plausible that it must be some sort of 
physical obstruction to ionic motion that gives rise to a slower effective diffusivity for Na"*". 
We then considered that physical obstructions to ionic motion can have the overall effect 
of increasing the time spent by each Na"*" in the nanospace thereby increasing the value of 
[Na''']ns- Ever since first reported measurement of Na"*" diffusivity in muscle cells[l8], it was 
hinted that it is physical rather than chemical interactions that cause a retardation of the 
ionic (and non-ionic for that matter — sorbitol and sucrose specifically) intracellular diffusion. 
Our ultrastructural observations of nanospace spanning electron opaque structures support 
this idea (Fig. HI) and agree fully with earlier reports of similar structures in SMC [21]. The 
introduction of such physical obstruction in the simulations seems then more than warranted 
and very plausible. Moreover, if the "tortuosity factor" , introduced by Kushmerick, arising 
from physical interaction is indeed important in the bulk cytosol, then it would be even 
more so in the restricted PM-SR nanospaces. While the introduction of physical obstacles to 
ionic motion in our simulations does produce higher [Na''']ns gains (Fig. [6] and [7]), to obtain 
a quantitative agreement with the observed level of [Na"*"] during the transients we need to 
hypothesize further that these junction spanning structures are not distributed at random, 
but rather form an organized barrier around a given Na"*" source within the nanospace. Our 
simulation results in this case (Fig. [9]) agree quantitatively with the observed [Na"*"] tran- 
sients, giving a Na"^ of the same order of magnitude as the experimental measurements. The 
left panel plot in Fig. [9] also shows that the [Na"*"] with this configuration of pillars will 
reach steady state with a time constant that is at least ten times longer than in the case of 
randomly distributed pillars. This is another feature supporting our hypothesis that these 
transients are formed in the PM-SR nanospaces with the aid of a relatively tight fence of 
pillars around the Na"*" source. In fact, the measured transients are reported to be very long 
lived with a ramp up time to steady state of about 30 s, a time scale that the simpler versions 
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Figure 8: To-scale x,?/-projection of model nanospace including 14 10-nm-radius cylindrical 
pillars positioned in a 5-nm-radius ring around a Na"*" source in the centre, thereby giving a 
1% porosity for ion passage. 

of the model could not reproduce at all. 

It therefore appears that two features are essential if these transients must occur within 
PM-SR nanospaces: there must be physical obstructions to Na^ motion, forming an orga- 
nized barrier around each Na^ source (TRPC6 channels) and NCX must be localized near 
a TRPC6 within such barrier to be able to sense the high [Na^], reverse and allow Ca^"*" 
entry. The second of these two features is suggested by a number of observations. There are 
studies supporting the idea that TRPC channels and Na'''/Ca^''' exchangers are functionally 
and physically linked and form an important signalplex in several different systems [29]. Still 
other studies indicate that at least certain kinds of TRP channels are found in close prox- 
imity of caveolin[30l [3T] . Combined with our earlier observations that NCX tend to crowd 
near the necks of caveolae|14j. this makes a strong case for a physical association between 
TRP channels and NCX. The first of the above mentioned desirable features — patterns of 
physical obstructions within the nanodomains — is a harder matter to study experimentally 
and one of the current /future directions our laboratory is following. 

To further study the issue of the localized [Na"*"] transient duration requires more com- 
putationally demanding simulations, which we are currently tackling. We can however make 
a qualitative argument to suggest how the time scale of the observed transients can also be 
supported by our model. We have so far only implemented one Na"*" source, or one TRPC6 
channel. However, we have observed earlier that during activation of rabbit inferior vena 
cava SMC there could be about 16 NCX in the junctions [13]. If we conjecture that TRPC6 
and NCX must be in physical proximity, as we speculated above, it is consistent to assume 
that there may be more than one TRPC6 in each junction functioning as a Na^ entry gate, 
and perhaps even as many as there are NCX. 

Other transporters such as NCX and NaJ^/K^ ATPases (NKA) certainly also play a 
role in shaping the PM-SR nanodomain [Na"*"] profiles that drive NCX-mediated Ca^"*" entry 
during VSMC activation. In the simple model we presented in this article, we focussed 
on the role of TRPC6 since it is by far the larger capacity transporter of the TRPC6, 
NCX, NKA trio (millions of ions per second, hundreds per second and tens per second, 
respectively). This first approximations quantitative model supports the idea that localized 
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Figure 9: [Na''']ns within a model nanospace containing 14 pillars "strategically" positioned 
around a Na"*" source at the centre and making a 1% porous fence, as in Fig [HI Left panel: 
computed [Na^]ns as a function of time. Right panel: concentration profile within nanospace 
at t = 0.0003 s. The red and the grey curves and dots refer to data displayed in Fig. [3] and 

[Na"*"] elevation transients take place at PM-SR nanodomains and strongly suggests that 
physical impediment to ionic mobility of Na"*" is also a necessary factor for their generation. 
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